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Abstract: We present an accurate and efficient electromagnetic particle-in-cell (EMPIC) 
algorithm on unstructured grids for the analysis and design of axisymmetric slow-wave structures. 
The use of unstructured grids allows for more fidelity in the modeling of micromachined 
geometries. The use of a reduced dimensionality algorithm decreases the computational costs 
significantly and enables its integration as a forward engine into a design loop. Special gather and 
scatter steps are employed to yield accurate beam plasma dynamics and exact charge conservation 
on unstructured grids. We provide numerical examples involving travelling-wave tube amplifiers 
designed to harness bunching effects arisen from Cherenkov radiation from plasma electron 
beams. 
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** Introduction 


“* Particle-In-Cell (PIC) algorithm for plasma simulations: 
Field solver, gather, relativistic particle pusher, and 
scatter 


“* Modeling of beam-field interactions in vacuum electronics 
devices (VED) 


% Numerical Examples (TWT amplifier and BWO) 


% Concluding Remarks 


ee Electromagnetic Particle-in-Cell (EM-PIC) Algorithms 





% Vlasov equation i Maxwell-Vlasov S tem 
= df, (r(t),v(t),t)/dt = 0 (collisionless) i D 


= Newton'sLaw 





O fp > 2 dp = > = > eo 
2 DET mo |[E+DxB]-V5f, =0 — OfMotion 
= Shape factor fp e presens spatial AQUA W A 
for an ensemble of particles (number density Equations 


of particles), described by a smaller set of “super- 
particles” (coarse graining of the phase-space). 
4 steps for each time-update 
% Electromagnetic particle-in-cell (€EM-PIC) algorithm 
= Numerical approach to solve the MV system 
= Four major steps at each time-update ————> 


f i ela = D de te 


“* Applications based on EM-PIC simulations 
= Plasma physics and astrophysics 
= Vacuum electron devices (VED) 





Particle- 
Update 
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+ Vacuum Electronics Devices (VED) 

= VED are essential for high-power microwave (HPM) sources in variety of applications 
including pulsed radar systems, RF amplification, and others. 

= Slow-wave structures (SWS) are often used to slow down the phase velocity of waves and 
excite Cerenkov radiation. 

=" Most VEDs are cylindrically axisymmetric (invariant along œ). 
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% Micromachined Slow-Wave Structures (SWS) 

= Micromachining desirable for high performance operation at microwaves and, especially, 
terahertz frequencies. 

= Either structured grids with conformal capabilities or unstructured grids (general solution) 
should be used to correctly capture the geometry of micromachined SWSs. 

= Long-standing challenge: Conventional EM-PIC algorithms on unstructured grids do not 
attain charge-conservation. 
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% Main contribution of this work 
= We introduce a EM-PIC algorithm for circularly symmetric VED which is (1) based on 
unstructured grids and (2) attains charge conservation from first principles. 
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“* 3D cylindrical axisymmetry VED 


= (1) Problem geometry, (2) EM fields, and (3) sources are cylindrically axisymmetric os =0). 


p 






ai 00, = R(z) (PEC) 


N 


vacuum Ep, Ez, Bg inQ — 


z J 


from charged particles 
to charged rings 
(invariant w.r.t. d) 


% Two useful constraints to reduce the dimensionality of 3D cylindrically axisymmetric 
problems 
1. Reduction of the exterior derivative d:d = dp + dz 


2. Consideration of only TM eigenmodes with m = 0 (m is an index for azimuthal modes) 


id Reduction of Problem Dimensionality 





% Differential forms expressions % Hodge star operator for fields and sources* 
for dynamic variables with a non H =x B J —* els 
orthonormal basis set q IS =i døndz* hdp ndh 
== = * Z * 
ide dedi LHS = Hodo p zP 
RHS =k a (Bydz A dp) RHS =x (J. dp + J,,dz) 
= -1p d d = yt B d 
D =x; E 
B = Bgdz ^ dp LHS = D,do ^ dz + D¿dp ^ dọ CERN 
H = Hydd RHS =*, (E,dp + E,dz) = a = pe j dp A dz 
= en lE E,(d d —* Pos) = 


= €(p)(E,dd Adz+E,dpAd 
J.=).,dp +),,d2 MS 40 nc + Bet ada) 


O, = Pv, 
J =] do Adz +J;dp ^do 
Q = p ,dp Ado A dz 

= We introduce an (1) artificial inhomogeneous medium with dependency on p and (2) rescaled sources. 





* Vey much akin to “transformation optics”: Teixeira and Chew, J. Math Phys. 40 (1), 169-187 (1999). 6 
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A Reduction of Problem Dimensionality (cont’) 





“* 3D axisymmetric problems w.r.t. geometry, fields, and sources can be modelled by 2D equivalent 
problems by TE®-polarized fields, artificial inhomogeneous medium, and scaled sources. 


p «* The continuous ring of particles 
along @-axis matches a particle 
at the meridian (poz) plane (arbitrary œ) at the meridian plane. 











Z 
“ring” of 
Z particles 
fields are TE?- artificial scaled sources P a particle at the 
polarized inhomogeneous AE 
. L 
medium 


A hbslis Mixed € — B FETD Scheme for (2+1) Theory 





de=-% dH="4+J, dD=Q, dB=0 


y spatial discretization by expanding variables as 


In primal mesh In dual mesh 
N ~ = 
E = Dj Erw (1—form) H= A H, 246 (0 — form) 


1 
N2 2 
B=}; 4B, °w? (2-—form) D= YN DP" (1-— form) 


at = p 
J.=YE 23, wi (1- form) JI=Y2,J, 07 (1- form) 


Q, = E, Qur wr (0 — form) Q= y Tee (0 — form) 


l=1 l l 


y Paring, generalized Stokes Theorem, and leapfrog time scheme 
Mixed € — B FETD Scheme: Kim and Teixeira, IEEE T-AP, vol. 59, p. 2350 (2011) 


(B]"*2 = [B]"-2 + At(Deun] « [EI 





[ke]: [E+ = [xe] -[E]" + At (Beun ea [1] p1? ) 


\ | 
A J : 
VA AN 


Q 
Z 
K : primal mesh 


K : dual mesh 


Discrete Hodge matrices 
[kelij = fe Ewi n w} 
— — 2 2 

[*,-1 ij Je y (pjwf A+ “i 
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Rao ven l 
e Gather and Particle Pusher 





% Gather 
= EM field values are interpolated at the instant 
positions of particles by using interpolants of 
Whitney forms 


© Ep = Lia EWE) 










1 1 1 1 n face-wise 
oNts >n-5 n+> n>>5 n IE 
a Bnp p _ y B; +B, | ye (22) E2 A constant 
p 2 ,=1 2 l p 
% Relativistic particle pusher ES 
A 

= Lorentz force and Newton’s law of motion t Ep © Bp 
= Updates of particles’ position and velocity | 


Relativistic factor, y = 1/ 1 (|? lic)” srt? sr 
= š — — p = p B qp =m >i = 
p Y ap = Tinga LED + Bp X Bp 


= Boris or Vay algorithms 
= Beam focusing system (BFS) Yy DP =p 2 
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% BFS prevents particles to radially deviate from electron beams due to the repulsive forces 


= It computes effects by self-fields and BFS and neglects @ motions due to continuous rings of particles 


ea neglecting the motion of @ due 
to continuous ring of particles 


n+1 


T p! 
4 nl 


N| e 














Pp, Zp 
le a 
“TNL NN VY 
E 0.6 - N i 
SRA the particle trajectory 
0.2 - at the meridian plane 
O' : i : t 
0 0.5 1 1.5 2 
z [m] 





(Projection of the 3D particles’ 
(3D particles’ surface trajectory under BFS) trajectory into the meridian plane) 10 
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A Exact Charge-Conserving Scatter 


Discrete continuity equation 





discrete charge 1] -[Q,]"*2 — [7] -[Q,]” 1 discrete current 
Qin; = q(x, pwr.) At Dee, = Ay Ab pwe,) 





TTA + p2 + + p2 + + Ps + + ps + 
2 n E @ - Ps af _ Ps E P3 4) E (z vite) = _ $ pe 2 Pa (asas — 2028) — ATA agas — 228) | 
Mt £ [oral — PsAq + (02 + pa (Al -11)] ~The [ppal — psAz + (pz + py — Ai)| 
Rate of charge variation at n4 Net currents flowing outward from n, 
n+1 


*x,N3 n+1 — n 


x,M4 x,N4 





difference 


— 





No 


osas 


“areas weighted by p” 


Example 1: Vacuum Diode 


L} Proposed charge-conserving (CC) scheme vs. conventional non charge-conserving (NCC) scheme 
Anode 


Particle dist. 
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Example 2: Traveling Wave Tube (TWT) Amp. 







Electron beam, input RF signal, SWS 
Sinusoidal Corrugated Circular Waveguide (SCCW) i | y 
Ry = 0.04, €g = 0.025, p¿ = 0.02 SOE 
Electron bunching due to forward Cerenkov radiation 2 ea 
Confinement of the electron beam by BFS 






Band gap 


D 


Frequency [GHZ] 
O1 


TM, modei A N 


coupling point 


3 
4 0.5 1 15 2 
B, P/T 


Dispersion relations for SCCW Distribution of the electron beam 13 


A ES Example 2: Traveling Wave Tube (TWT) Amp. (cont?) 
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% Residual analysis for discrete Gauss law (DGL) = Normalized residuals for DGL 




















IH n — n 
DGL: IDaiv | f [*e] [E] — [Q] [Q] 
NR5cL, = 1- STA SM TEM) 
= The residuals for DGL at any two successive time-steps 2 ¡cal aiv) ¡(Za l%el xl lk) 
remain the same, in other words, if initial conditions have | | 
lDaiv| -[x-] - [E]? = [Q,]9, DGL for all time-steps is verified. 10° - t = 400,000At i 
Premultiplying |[Dgiv| - [ke] into DAL, using [Daiy| - |Dcun| = 0, and 
rearranging the Bid ac E 
10° - 
a [E]"**—[E]” n+> — e. 
[Daiv| [el ( At )= o =[Daiv) l Da] 2. ` D residual level 
Substituting DCE into above equation and rearranging it give z / 
= -10 
E 10 F a - 
[Daiv] - [x] - [E]"*? — [Q,]"*1 = = |Daiv| - [ - [x] - [E]” — [Q,]”. double precision floor 
By induction, 
[Daiv] + [xe] - [E]" — QJ” = [Baiv]: e]: [E — [Q.]° 0 E IN 
— oo, a — o 200 400 600 800 1000 


n 0 ; : 
res res Nodal index, / 


Example 3: Backward Wave Oscillator (BWO) 


= Small perturbation by the electron beam generates RF signals in the positive feedback system (BWO). 
= Rectangular corrugated cylindrical waveguide : Rin = 0.04 [m], Rout = 0.045 [m], pe = 0.045 [m] 


" V, = 500 [kV] yields relativistic electron beams with v, = 0.83c. 
= Electron bunching due to backward Cerenkov radiation 
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= New electromagnetic particle-in-cell (EM-PIC) algorithm developed to solve 
Maxwell-Vlasov system on unstructured grids. 

= Exterior calculus of differential formst was used to derive an exact charge- 
conserving gather and scatter schemes on unstructured grids?. 

= SPAI algorithm was used to accelerate the field update”. 

= Relativistic particle pusher has been implemented by using Boris algorithm. 

= Successfully applied to simulate plasma-beam dynamics in axisymmetric 
VED devices, including travelling wave tube (TWT) amplifiers and backward 
wave oscillators (BWO). 
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